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The spin degree of freedom in solids offers opportunities beyond charge-based elec- 
tronics and is actively investigated for both spintronics and quantum computation. 
However, the interplay of these spins with their native environment can give rise to 
detrimental effects such as spin relaxation and decoherence. Here, we use an all- 
electrical, lateral GaAs spin valve to manipulate and investigate the inherent nuclear 
spin system. Hanle satellites are used to determine the nuclear spin relaxation rates 
for the previously unexplored temperature range down to 100 mK, giving Ti times 
as long as 3 hours. Despite metallic temperature dependence of resistivity, the ob- 
served relaxation rates show a sub-linear temperature dependence. This contrasts 
the Korringa relaxation mechanism observed in metals but is not inconsistent with 
hyperfine-mediated relaxation in a disordered, interacting conductor not far from the 
metal-insulator transition. We discuss possible relaxation mechanisms and further 
investigate inhomogeneities in the nuclear spin polarization. 
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Spintronics is a promising vision for next generation 
electronics initiated by the spin transistor proposed 
by Datta and Das @, HI- Spin-based phenomena have 
become building blocks for both logic and memory de- 
vices 0, HI- Furthermore, the proposal of an electron- 
spin quantum computer [6] has opened up a very active 
field of research aimed at understanding and controlling 
spins in semiconductor quantum dots where the nu- 
clear spin bath of the host material has been identified 
to be of major importance. The interplay of the electron 
and nuclear spin systems in GaAs nanostructures leads 
to many interesting effects, including dynamic nuclear 
polarization (DNP) intricate mutual feedback 

Ojj], electron spin decoherence and further allows 

coherent manipulation of the nuclear system Weak 
coupling leads to long nuclear relaxation and coherence 
times which makes the nuclear spins themselves interest- 
ing for logic and memory applications. 

Here, we use a lateral spin valve device [lU [l6[ as a 
model system to investigate the GaAs nuclear spin envi- 
ronment. A spin valve offers an efficient way to i nject , 
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manipulate and detect electron spins in GaAs 
Electron sp in p olarization can be transferred by a DNP 
process (2fjl-[23j to the nuclear system, mediated by the 
hyperfine contact interaction via spin flip-flop processes. 
The spin valve thus allows us to manipulate and investi- 
gate the native nuclear spin environment by all-electrical 
means, rather than using optics 24J, making easily ac- 



cessible the low-temperature regime T <C 1 K which was 
previously not explored. The nuclear spin polarization 
acts back on the electron spins as an effective Overhauser 
field -Bat [HI causing electron spin precession. This re- 
sults in clear features in spin transport measurements, 



including a depolarization peak [rjj, l!8|, |26( and satellite epilayer charge carrier density of n = 5 x 10 cm 



field cancels the internal Overhauser field Bjv 0. The 
depolarization peak saturates already for relatively small 
nuclear polarization, but the Hanle satellites can easily 
be used as a sensitive probe of Bjv, in particular for large 
polarizations and down to the lowest temperatures. 

Employing a novel, all-electrical pump-probe cycle, 
we measure the nuclear spin relaxation time using the 
Hanle satellites over a broad temperature range from 
100 mK to 20 K. At high temperatures T > 1 K, the mea- 
sured Ti times are in good agreement with NMR data 
in similar systems 27J, |28[ . At the lowest temperatures 
T ~ 100 mK, we observe T\ times as long as 3 hours. 
Interestingly, despite metallic resistivity of the spin valve 
devices, we find that the temperature dependence of the 
nuclear spin relaxation rate follows a power law T o e 
rather than the Korringa law T^ 1 oc T |[20j ] naively ex- 
pected for a metallic system. However, we argue that 
the combination of disorder and electron-electron interac- 
tions (3^1 in a degenerate semiconductor not far from the 
metal-insulator transition can lead to a deviation from 
the standard Korringa law within the hyperfine-mediated 
nuclear spin relaxation mechanism applicable here. Fur- 
ther, the effects of inhomogeneities in the nuclear spin 
ensemble are investigated. 

The spin valve device used in this study consists of five 

6 nm thick Fe bars on a c(4x4) reconstructed surface of a 1 

micron thick epilayer of n-doped GaAs, see Fig.[Tp. The 

three contacts in the center have widths of 6, 2 and 1 /jm 

and edge-to-edge distances of 3 and 4.5 /*m, respectively. 

A 15 nm thick highly doped surface layer allows efficient 

spin injection and detection across the Schottky barrier 

formed at the metal-semiconductor interface. Within the 

subsequent 15 nm, the Si donor density is reduced to the 
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peaks in the Hanle geometry when the applied magnetic Further device preparation details are described in |18j . 
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FIG. 1: Signatures of nuclear spin polarization in spin 
valve and Hanle measurements, a, Spin valve measure- 
ments at 4 K showing the non-local voltage Vnl versus Bx 
applied along the Fe bars. An offset Vb is subtracted. Red 
curve is negative sweep direction, blue positive. Each pair 
of bold arrows designates the detector and injector magne- 
tization direction, respectively, b, Device and measurement 
schematic, c and d, Hanle measurements at 4 K. Displayed is 
Vnl (in d as color scale) versus Bz applied perpendicular to 
the sample plane for various Bx- The same parabolic back- 
ground is subtracted for all Bz sweeps (0.34 mT/s). Dashed 
lines correspond to cuts shown in c. We use the satellite peaks 
seen for Bx < as a measure of Bn, see text. 



The measurement scheme used drives currents from 
the injector 2 to the 100 /im distant contact 1 while 
a non-local voltage Vnl is measured between contacts 
3 and 5 outside the charge current path (see Fig.[Tp). 
Vnl is detected via a lock-in technique using a small ac- 
modulation on top of a dc- injection current. Electron 
spin polarization is efficiently injected into the semicon- 
ductor 17, 3l|, where it diffuses away from the place of 
injection. The electron spin diffusion length Id is several 
microns, larger than the distance to the detector contact 
[i~8| . The measurements are performed in a dilution re- 
frigerator equipped with a home-built 3-axis vector mag- 
net, allowing us to determine the magnetization direction 
of the iron bars to better than 1° by rotating the mag- 
netic field during continued spin valve measurements. 

Spin valve measurements are shown in Fig.fT^. Clearly 
distinguishable voltage levels for parallel and anti-parallel 
alignment of injector and detector magnetization are ap- 
parent. Further, the non-local spin signal appears sup- 
pressed around Bx ~ [l7j . as the electron spins start 
to precess and depolarize due to Bn ^ [l8j |. When Bn 
is strong enough to completely wash out the polarization 
of electron spins on their diffusive path to the detector, 
the depolarization dip saturates and becomes insensitive 
for larger Bn- Therefore, the depolarization dip cannot 
easily serve as a sensitive detector for large nuclear spin 
polarizations. 



To overcome this limitation, we turn to Hanle mea- 
surements [TBI, , where an additional external field Bz 
is scanned perpendicular to the sample plane, see Fig. 
[TJ; and d. The injected electron spins precess, diffuse and 
dephase in the perpendicular field Bz- Around Bz ~ 0, 
the spins do not precess much and reach the detector 
mostly unchanged, giving a peak in the spin signal (Hanle 
line-shape). For the following, both injector and detec- 
tor were initialized along the positive x-direction and are 
kept parallel throughout. For Bx < 0, two satellite 
peaks appear (Hanle satellites) [18| as seen in Fig. [TJ; 
and d, displaying a (partial) recovery of the spin signal 
and hence a suppression of spin precession and dephas- 
ing [32|, [33| . The satellite peaks thus indicate that Bz is 
effectively canceled by the z-component of Bn- 

Theory estimates the Overhauser field B^r in steady- 
state and for large B ext — \H e xt\ 3> Bl, where Bl 
1 mT is the rms field seen by a nuclear spin, 



as 
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where b° N is the maximum nuclear field obtainable with 
DNP and \Sx\ < 1 is the resident electron spin polar- 
ization along the x-direction. Following eq. ([T]), Bjv is 
either parallel or antiparallel to B ea;t and in particular, 
B^ points opposite B ex t for Bx < and S > 0. Thus, 
the external field vector 'Bext can cancel the Overhauser 
vector Bj\r. This leads to a recovery of the electron spin 
polarization and appearance of the satellite peaks, which 
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can therefore be used to measure B^r 
Bz ^> Bx,By, we have Bn ~ Bz on the satellite peak. 
Nuclear fields achieved here are ~ 50 mT, about one per- 
cent of the ~ 5 T for fully polarized nuclei in GaAs. 

The steady-state situation assumed in eq. (JXJ) is not 
typically reached in our measurements, since the ramp 
rates employed are fast compared to the time scales of 
the nuclear spins. Perfor ming a time average (-)t of the 
magnitude of eq. (JXJ) gives [18| 
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Note that the largest contributions to (-) t arise around 
Bz ~ and therefore B ext ~ Bx, where the equilibrium 
Bn = b N Sx is maximal. The resulting average nuclear 
field Bn and therefore satellite peak splitting is linear in 
Bx, as seen in Fig.QJl [18]. Further, significant broaden- 
ing of all three peaks is visible for increasingly negative 
Bx- For the following, we fix Bx — — 1.5mT as a com- 
promise between peak splitting and broadening. 

Nuclear spin relaxation times T\ are measured using 
a pump, decay and probe cycle as sketched in Fig. [5^. 
First, the nuclear polarization is built up by injecting a 
large current for efficient DNP while continuously sweep- 
ing Bz back and forth ('initialize'). Nuclear spin polar- 
ization is gradually accumulated over time, see Fig.[5p, 
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FIG. 2: Nuclear spin relaxation determined from 
Hanle measurements using a pump, decay and probe 
cycle, a, Measurement schematic, b, Initialization: Hanle 
Bz sweeps (0.3mT/s) as a function of time with Ida = 20 fiA. 
Satellite peaks spacing grows and eventually saturates after 
typically one hour. Trace begins at Bz = +75 mT, sweeps 
to Bz = —75 mT and then continues to alternate sweep di- 
rection, c, Hanle Bz probe-traces (0.9mT/s) as a function 
of delay time r, here starting from a slightly more polarized 
nuclear state compared to the end of b. Decay of the satel- 
lite peaks is clearly visible. A parabolic background identi- 
cal for all delays r was subtracted, d and e, Log-plot of 
Overhauser field Bn (crosses) - extracted from satellite peak 
positions such as in c - as a function of r measured at 4.2 K 
in d and 170 mK in e. Blue data is from negative- -Bz, red 
from positive- Bz satellite peaks, with single-exponential fits 
(solid lines) giving excellent agreement. Very long Ti times 
are obtained, as indicated. 



eventually reaching a steady state, typically after an 
hour. As will be shown later, the continuous sweeping 
has the effect of homogenizing the nuclear spin system, 
producing sharper satellites. The asymmetry, alternating 
positions and alternating widths of the satellite peaks is 
a consequence of the time average in eq. @ , arising from 
the alternating sweep directions and efficient DNP when 
Bz ~ followed by some decay and homogenization for 
\Bz\ 3> 0. After initialization, the nuclear polarization 
is allowed to decay for a time r with Idcac — and Bz 
ramped to zero ('decay', keeping Bx = — 1.5 mT fixed). 
Subsequently, a fast Hanle scan is performed with only a 
small ac-current (no dc-current) applied to avoid further 
DNP ('probe'). Repeating this cycle for various delays r, 
data sets reflecting the decay of the nuclear Overhauser 
field Bn over time are obtained, as shown in Fig.(2£. Note 
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FIG. 3: Temperature dependence of the nuclear spin 
relaxation rate. T 1 _1 (T) is shown between 20 mK and 
20 K measured in two consecutive cool downs (open and solid 
squares) on the same sample, all for Bx = — 1.5 mT. Er- 
ror bars are estimated from repeated measurements under 
identical conditions. The solid (red) line is a power law fit 
1/Ti oc T a with a = 0.6 ± 0.04 (excluding the data below 
100 mK from the fit). As a comparison, the slope of the Ko- 
rringa relaxation law 1/Ti oc T applicable in metals is added 
(dashed line). Upper inset: Bx dependence of the nuclear 
Ti-rate at T=10 K with theory fit (dashed curve) to eq. |3|. 
see text. Lower inset: temperature dependence of the resis- 
tivity from van der Pauw measurements on the same GaAs 
wafer, indicating metallic behavior for T < 10 K. 



that the overall measurement time is 10-30 hours due to 
the very long nuclear time scales and repeated initializa- 
tion for each r. 

By fitting Lorentzians to the satellite peaks, we obtain 
Bn (peak position) as a function of r, as shown (crosses) 
in Fig.di at T = 4.2 K and Fig.d? at T = 170 mK, 
for both positive (red) and negative (blue) Bz satellites. 
The small difference between the two satellite positions 
is again a result of eq. @ and the sweep direction. From 
exponential fits we get excellent agreement with the data 
and nearly identical T\ times for the two satellites. We 
note that Bn(t) clearly follows a single-exponential de- 
cay over the measurable range oi Bn (about one order of 
magnitude in Bn) at all temperatures. In addition to the 
decay of Bn we observe sharpening of the satellites with 
increasing delay times, indicating increasing homogene- 
ity of the nuclear spins with time. At high temperatures 
> 1 K, the T\ times obtained here are in good agreement 
with previous T\ measurements by NMR on GaAs with 
similar doping (27ll2^|. 

The temperature dependence of the nuclear spin re- 
laxation rate is shown in Fig.FJ] on a log-log plot for 
two cool downs (open and closed squares) of the same 
sample. Measurements of a second sample (not shown) 
fabricated on the same wafer give very similar results. 
Both ac and dc currents were experimentally chosen to 
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avoid self-heating. However, in the refrigerator we used, 
sample temperatures saturate around 100 mK due to in- 
sufficient thermalization, causing the relaxation rates to 
saturate at 100 mK. Nevertheless, at the lowest temper- 
atures, very long T\ times exceeding 3 hours are found. 
Motivated by the approximately linear temperature de- 
pendence in the log-log plot ranging over two orders of 
magnitude in T, we fit a power law 1/T\ oc T a for 
100 mK< T < 10 K and find a = 0.6 ± 0.04. The data 
points at T > 10 K are excluded from the fit due to well 
known phonon contributions above 10 K [27j . 

In simple metals, the Korringa law 1/Ti oc T is ex- 
pected |29|. Indeed, this is observed in much higher 
doped bulk GaAs (n = 2xl0 18 cm~ 3 ) j27| measured with 



NMR above 1 K. These T x data from Ref. [27j are added 
to Fig. [3] for comparison (dashed line, extrapolated below 
IK). This clearly contrasts with our results, which show 
a sub-linear temperature dependence extending well into 
the previously unexplored temperature range below 1 K. 

To shed more light on the unexpected temperature 
dependence of T-j - , we first investigate the resistivity 
p. The lower inset of Fig. [3] shows p(T) from van der 
Pauw measurements done on separate samples from the 
same wafer. A weak T dependence with metallic be- 
havior (dp/dT > 0) is seen for T < 10 K, as expected 
for the present doping of 5 x 10 16 cm~ 3 , above the well- 
known metal-insulator transition (MIT) in GaAs at ~ 



2 x 10 16 cm- 3 [34j,|35|. The charge density at 4 K is the 



same as at base temperature within measurement error. 
Further, a perpendicular magnetic field at low T has no 
significant influence on the resistivity. Therefore, though 
the density is only a factor 2.5 above the MIT, the re- 
sistivity data shows clear metallic behavior, lacking any 
hints of incipient localization. In addition, we have con- 
firmed that the highly doped surface layer does not sig- 
nificantly contribute to transport apart from facilitating 
the Fe contacts. The interaction parameter r$ — Ec/Ep 
is about 0.6, with Fermi energy Ep = 7.4 meV and aver- 
age Coulomb energy Ec — 4.1 meV, indicating that the 
samples are approaching the interacting regime r$ ^ 1. 
Further, disorder is quite strong: kp£ ~ 1.7, with Fermi 
wavelength Xp = 55 nm and transport mean free path 
£ = 15 nm for T < 10 K. Therefore, the epilayer be- 
haves like a degenerately doped semiconductor showing 
clear metallic behavior, in the somewhat interacting and 
strongly disordered regime. 

We now discuss the temperature dependence of the nu- 
clear spin relaxation rate. First, we exclude phonon con- 
tributions since these have been shown to be relevant only 
well above 10 K and would result in a quadratic temper- 
ature dependence T^ 1 oc T 2 27, 3(J. Also, nuclear spin 
relaxation by paramagnetic impurities is know to be very 
weak in doped GaAs [27[. Next, we consider nuclear spin 
diffusion out of the 1 pm thick epilayer. This process is 
in principle temperature independent in the regime ap- 
plicable here, but combined with some other relaxation 



channel - e.g. T x 1 oc T a with a > - it could be- 
come dominant at low T. However, the nuclear spin dif- 
fusion length £ N = yJD N T x with D N = 10- 13 cm 2 /s ^ 
is smaller than the epilayer thickness, even at low T. Fur- 
thermore, this random- walk contribution is inconsistent 
with the clear single-exponential decay of Bn(t) which 
we find for all temperatures, also making it unlikely that 
the observed low-T saturation of T± is caused by nu- 
clear spin diffusion. Therefore, we estimate that nuclear 
spin diffusion is negligible here. 

Next, we consider the hyperfine contact interaction as 
a possible nuclear spin relaxation mechanism. In non- 
degenerate semiconductors, the Fermi energy is well be- 
low the conduction band edge, the mobile charge carriers 
follow a Boltzmann distribution, and the nuclear spin 
relaxation rate is Tf 1 oc VT 38], not far from the mea- 
sured T^ 1 oc T 6 . However, since the measured resistiv- 
ity p(T) does not display a thermally activated behavior 
expected for a non-degenerate semiconductor, this mech- 
anism is most likely not applicable for a homogeneous 
system. In presence of spatial inhomogeneities with non- 
degenerate (low density) and degenerate (high density) 
electron regions, a combination of Korringa (T^ 1 oc T) 
and Boltzmann relaxation (T£~ oc VT) could in principle 
lead to the observed T^ 1 oc T 6 . However, for the low- 
density regions to remain in the Boltzmann regime down 
to ~ 100 mK, the density would have to be several orders 
of magnitude lower than bulk, making this unlikely. 

To learn more about the mechanism of nuclear spin 
relaxation, we investigate the Bx dependence of J\ , 
shown in the upper inset of Fig. [3] at 10 K. Note that 
Bz = during the decay cycle of the Ti measurement. 
A clear reduction of relaxation rates is seen for increas- 
ing Bx, as expected for applied fields comparable with 
Bl , which is the local rms field acting on each individual 
nuclear spin, including nuclear dipole-dipole fields and 
electronic Knight fields. The theoretically expected rate 

is &Si 



r -i (m _ B* + 6(5/3)Bl 



B 2 + (5/3)51 



(3) 



with large-field rate a = T 1 1 {B > B L ). Note that the 
zero- field rate T^ 1 (B = 0) = 8a and the correlation pa- 
rameter 8 is ranging from 2 for uncorrelated to 3 for fully 
spatially correlated fields B^. Independent measure- 
ments give a very small B-field offset < 0.1 mT, which 
we assume to be zero here. We perform a fit constrain- 
ing 2 < 8 < 3 and obtain 8 = 3.Q ± 0.3, B L = l± 0.2 mT 
and a — (9.6 ± 1.5) x 10 _4 s _1 . Further, the dashed the- 
ory curve shows a best fit with 8 = 2 held fixed, which is 
clearly inconsistent with the present data. Therefore, Bl 
is spatially highly correlated (8 = 3) with local field Bl 
larger than estimated nuclear dipole-dipole fields alone. 
This suggests an electronically induced hyperfine Fermi 
contact interaction mechanism causing nuclear spin re- 
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laxation, due to extended electrons on a length scale 
much larger than the lattice constant. 

For metals or degenerate semiconductors, the hyper- 
fine contact interaction results in a nuclear spin relax- 
ation rate in the free-electron approximation given by 
[38| (Korringa mechanism) 



a 0.50 



1 



256tt 3 7 ^ 

9h 7e 2 



n 2 \m\ A X 2 -k B T , 
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with gyromagnetic ratio 7„ of the nuclei and j e of the 
electrons, electron spin susceptibility \ an d electron den- 
sity 77. 1 0(0) | 2 at the nuclear site. This mechanism is 
consistent with the measured Bx dependence indicat- 
ing strongly correlated fields Bl and the observed single- 
exponential decays of Bm{t). Beyond the free-electron 
model, Tj -1 would need to be calculated properly in- 
cluding the combined effects of disorder and interactions 
M, El El not far from the MIT. In lack of a Tf 1 the- 
ory in this regime, naively, a corresponding renormalized 
electron spin susceptibility appears in eg. Q and can de- 
velop a temperature dependence (43l446j |. This was also 
seen in other semiconductors for 



> 77 c above but not 
far from the MIT 47J or the Stoner instability. Here, 
X oc T~P with j3 = 0.2 ± 0.02 would be required to result 
in T^ 1 oc T 6 as measured, assuming no other T depen- 
dencies in eq. (H]). Such electronic correlations and x(T) 
are potentially caused by disordered, interacting mobile 
electrons, RKKY like effects or spin orbit coupling. In 
addition, some fraction of localized electrons in the tails 
of the impurity band can contribute to nuclear spin re- 
laxation [271 ] and can further cause exchange correlations 
with delocalized electrons, such as the Kondo effect and 
hierarchical formation of spin singlets [27, 47, 48 1. The 
density dependence of j5 would be interesting to investi- 
gate, indeed, but is beyond the scope of this study. As the 
free-electron regime is approached for increasing densities 
77 ^> n c , one would expect f3 — > 0, as already reported for 
larger densities [27j |. 

Finally, we investigate nuclear spin inhomogeneities 
apparent in Hanle measurements. Fundamentally, the 
decay of the electron spin polarization over the space be- 
tween injection and detection contact can lead to spatial 
decay and inhomogeneity in the induced nuclear spin po- 
larization and corresponding satellite peak broadening. 
However, the electron spin diffusion length is larger here 
than the injector-detector distance [lij . Further, given 
sufficient time, nuclear spin relaxation (and potentially 
diffusion) can counteract polarization imbalances, mak- 
ing the nuclear inhomogeneities sensitive to the details 
of the DNP scheme such as sweep times and Bz field 
conditions. 

Since the strongest DNP occurs around Bz ~ (see 
eq. @), we fix Bz =0 during initialization (rather than 
sweeping it), followed by decay and probe cycles as for 
the T, measurements previously. As seen in Fig.0^ and 
b, this results in significantly broadened satellite peaks 
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FIG. 4: Effects of nuclear spin inhomogeneities. a and 

b, are Hanle probe measurements from a pump, decay and 
probe cycle as in Fig. [2] but here with initialization done at 
Bz = fixed rather than sweeping. Though the satellite 
peaks are obviously much broader (than e.g. Fig.[5|, indi- 
cating stronger nuclear inhomogeneities, very similar nuclear 
relaxation times are extracted, c and d, In slow Hanle mea- 
surements (0.25 mT/s), additional satellite peaks become vis- 
ible (arrows in c) for Bx > 0.8 mT, see text. 



compared to Fig.j^t, consistent with the absence of a nu- 
clear spin homogenization mechanism as discussed before 
in Fig.[2t>- Despite the additional broadening, however, 
the extracted nuclear spin relaxation times (not shown) 
are identical within the experimental error to the sweep- 
ing Bz initialization in Fig. [21 Further, clear sharpening 
of satellite peaks is seen in Fig. [4] with progressing delay 
r, consistent with our previous observations in Fig.[2jD 
and c. 

Additional satellite peaks are present in Hanle mea- 
surements, see Fig.0J;, clearly visible but smaller and 
more broadened than the main satellites. Note that dou- 
ble satellites appear whenever the Hanle satellites are 
well enough separated and sharp enough, irrespective of 
current direction and sign of Bx , here shown for Bx > 
and reversed current in Fig.@jL The additional satellites 
have also been observed on a second device and seem to 
suggest two distinct species of electron and/or nuclear po- 
larization regions between injector and detector. Drifting 
electrons will also feel the spin orbit interaction, which 
might alter the effective magnetic field they are exposed 
to. Further studies are needed to elucidate these addi- 
tional peaks. 

In summary, the presented measurements on GaAs 
spin valves show the importance of the nuclear spin envi- 
ronment, which can easily be manipulated all-electrically 
via DNP at temperatures below 20 K. Extremely long nu- 
clear spin relaxation times exceeding 3 hours are seen. 
The nuclear spin relaxation rate has a rather weak, 
power-law temperature dependence T^ 1 cx T 6 , devi- 
ating from the standard Korringa law T± oc T in met- 
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als, but not inconsistent with nuclear spin relaxation via 
the hyperfine contact interaction in the interacting, dis- 
ordered electron regime currently lacking theory. Inho- 
mogeneities in the nuclear spin polarization were inves- 
tigated, and appeared not to affect the spin relaxation 
rate. 
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